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Abstract 

Mesoscopic effects associated with wave propagation in spacetime with metric stochas- 
ticity are studied. We show that the scalar and spinor waves in a stochastic spacetime 
behave similarly to the electrons in a disordered system. Viewing this as the quantum 
transport problem, mesoscopic fluctuations in such a spacetime are discussed. The con- 
ductance and its fluctuations are expressed in terms of a nonlinear sigma model in the 
closed time path formalism. We show that the conductance fluctuations are universal, 
independent of the volume of the stochastic region and the amount of stochasticity. 
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1 Introduction 



1.1 Metric fluctuations in semiclassical gravity 

Since Einstein-Hilbert action is insensitive to fluctuations near the Planck scale, we expect 
large fluctuations in this regime to occur which may require the reconsideration of the concept 
of microstructure of spacetime itself. The concept of metric fluctuations thus introduced 
originally by Wheeler |jl|, |^ have been studied and modifled in various different contexts. 
Spacetime foams or other string theory motivated microstructure of spacetime can be treated 
as possible stochastic sources 0. Since those fluctuations become dominant only near the 
Planck scale and cannot be directly observable, it is important to identify their possible 
influence on the matter flelds propagating in such a background with the energy much lower 
than the Planck value which is possibly detectable by high energy astrophysical observation. 

Since quantum gravity is still an unsolved problem, in the conventional semiclassical 
gravity, a spacetime is left unquantized 0]. The recent study in semiclassical gravity, how- 
ever, reveals many pathological features in this approach. One of the most serious problems 
is the violation of positive energy theorem due to large quantum fluctuations which leads 
to the violation of causality by allowing the creation of traversable wormholes It is also 
noticed that the positivity can be recovered by imposing the additional smearing in the case 
of Minkowski spacetime ^ . This smearing, originated in the microscopic quantum fluctua- 
tions, may also cure other problems such as initial or blackhole singularity problems. Other 
studies also suggest that the small stochastic fluctuations of spacetime metric lying on the 
deterministic background spactime are not only the useful phenomenological modiflcation of 
semiclassical Einstein equation but also the inevitable consequence of the more fundamental 
processes or of the backreaction induced by matter fluctuations [0, In an astrophysical 
context, squeezed states, evolved from primordial gravitational waves by parametric ampli- 
flcation during the cosmological expansion [^], are considered to be one of the origins of the 
stochastic gravitational waves believed to exist in the present universe. Whether or not the 
trace of such primodial fluctuations as the quatum noise lies within the detectable range of 
the Laser Interferometric Gravitational Wave Observatory detector (LIGO)[0 is the curious 
question under debate [|ll|, |l^ . 

1.2 Semiclassical gravity and mesoscopic physics 

Semiclassical gravity, though far different in energy scale, shares many common features 
with mesoscopic physics |T3| , The quantum transport properties of metallic systems are 



known to be divided into several different regimes depending on the qualitatively different 
contribution of scatterers. For the length scale less than the mean free path Im, the wave 
propagates ballistically similar to the free coherent wave; for the length scale larger than 
the coherence length Lcoh, the scale at which the mutual coherence of waves is lost due to 
inelastic scattering, the classical Boltzmann transport theory is valid. In the mesoscopic 
scale Im < L < Lcoh, multiple scattering has to be taken into account and the coherence be- 
tween different paths becomes important. The effect of the environment has to be taken into 
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account and the dissipation and decoherence due to inelastic scattering play an important 
role. Furthermore, in the mesoscopic regime close to L^oh) the quantum-classical correspon- 
dence of the propagating wave becomes relevant; in the same regime close to /m, possible 
influence of the microscopic constituent of a medium will manifest itself on the transport 
properties of the wave. Owing to the recent progress in nanoscale technology [0, many 
phenomena in this regime are amenable to experiments. In light of this analogy, various 
effects associated with the electromagnetic wave propagation in the Friedmann-Robertson- 
Walker universe and the Schwarzschild spacetime with metric stochasticity were studied in 
T6| based on the formal equivalence of the Maxwell's equations in a stochastic spacetime 



with those in random media. Localization of photon and anomalous particle creation can 
occur in such stochastic spacetimes. 

In this paper, we show that the scalar and spinor fields propagating in a stochastic 
Minkowski spacetime can be considered as the electrons propagating in a disordered po- 
tential. The randomness couples to the frequency of the wave additively for the conformal 
metric fluctuations. Thus the effect of the stochasticity resembles that of a random potential. 
We use the closed time path formalism that allows us to study equilibrium and nonequilib- 



rium quantum field theory in the unified framework [|T^ [T^. The closed time path partition 
function is defined for the stochastic quantum system from which we obtain the effective 
interaction between fundamental fields in the form of a four point vertex after averaging 
over stochasticity. Diagonalization of the matrix Green function is employed. This is an es- 
sential step to make use of the nonlinear sigma model developed for the disordered transport 



problem |T^, ^ We use the Hubbard-Storatonovich transformation and write 

the collective excitations in terms of auxiliary fields. The conductance and its fluctuations 
are expressed by these auxiliary fields. 

The paper is organized as follows: In Sec. 2, we start from showing that the effects 
of conformal metric fluctuations on scalar and Dirac fields can be simply represented by 
the effects of a fluctuating mass. In Sec. 3, we develop the nonlinear sigma model to 
express higher order fluctuations of fundamental fields in terms of the correlation function 
of collective fields. The conductance associated with propagation of particles in such a 
spacetime can be defined analogous to the one in electric circuits by the Kubo formula as 
the correlation function of currents. Similar to the mesoscopic quantum transport problem, 
we show that the conductance fluctuations are universal, independent of the size of the 
stochastic region and the amount of stochasticity initially assumed. The amplitude of the 
conductance fluctuations is constant upto the leading order in the weak disorder expansion. 
In Sec. 4, the summary of results is given followed by brief discussions. 



2 Wave Propagation in Stochastic Spacetimes 
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2.1 Scalar and Dirac fields in stochastic spacetimes 

The Lagrangian density for the complex free scalar field in curved spacetime has the form 

C-s = ^r9[9,.d^<l>^d^<l> - {ml + iR)(f'l (1) 

where ^ is a dimensionless nonminimal coupling parameter. ^ = 1/6 is called conformal 
coupling and = is minimal coupling. In the presence of a slight amount of inhomogeneity 
in a flat spacetime background characterized by g^i, = rj^^ + h^^, the Lagrangian given above 
can be split into the flat space term and the perturbation term as £5- = Csq + Csf- 

Cso = ri^'d.cjy^dA - mlcj)^ (2) 

and Csi = ~Ts'^h^^, where Tg'^ is the stress-energy tensor. In a conformal coupling case 
= 1/6), the action Eq. (|lD is conformally invariant except for the mass term. Thus, a 
conformal type of metric fluctuations can be attributed to the effect of a fluctuating mass. 
Accordingly if we write the stress-energy tensor in the following form 

+ (f)^d^dy-^r]^y^n(f)-^T^'"'ml(j)^] (3) 

for an isotropic and conformal type of stochasticity (7^1/ = r]^ve"'^^^ ~ y]^u + ^ij.vv{x), where 
v{x) is a stochastic field, the interaction term can be simply written as 

Csi = -mlv{x)<l)\x). (4) 

In this case the total Hamiltonian is given by 

H = J d^x[7r^ + + ml(f)'^ + mlv{x)(f)% (5) 

or in a momentum representation, 

p <? 
Furthermore, if we restrict our case to elastic scattering, we have 

p ? 

where uj =p^ . This has the form equivalent to electron propagation in a disordered potential. 
The equation of motion in this model will be 

V'^(t)+{uj^ -ml)<P + vs{x)(t) = Q, (8) 
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where vs{x) = m|f (x). Comparing Eq. with the equation of motion for the electromag- 
netic wave |[I6[] 



V^(j) + uj^(f)-uj^vs{x)(f) = 0, (9) 

we see that the randomness couples additively to the frequency in the scalar wave, whereas 
it couples multiplicatively in the electromagnetic wave. Since the electromagnetic wave 
obeys Maxwell's equations which are conformally-invariant, the equation of motion is rather 
similar to that of wave propagation in random media. While for the scalar wave subjected to 
conformal metric fluctuations, the equation of motion has the form of massive particles in a 
random potential. We expect that the property of fields under the influence of more general 
types of metric fluctuations will possess both aspects. The difference of the form between 
Eq. (H) and Eq. (||) leads to the different density of states and transport properties. We 
also point out that although the restriction of the metric fluctuations to the conformal type 
simplifies the arguments significantly, many arguments in the rest of the paper is applicable 
to the more general class of fluctuations. 

The Lagrangian density for the Dirac field in curved spacetime has the form 



D 



1 _ _ ' 



(10) 



where ^7 = 7aV°, 7° = &"^7'^, and b = detfo'^^ for vierbein fields b^^. The Lagrangian given 
above splits into the flat space contribution and the perturbation term as C = Cdo + Cdi 
similar to the scalar field case: 

1 _ _ 
^Do = -^S^ - mDipilJ (11) 

and Cdi = —Tj^'^hf^,^, where Tj^'^ is the stress energy tensor that can be written as 

i r T 



rpflU 



^-ff'W^ -W^lj-f^^ . (12) 



For a massless Dirac field, the action Eq. ( pUj ) is conformally invariant. Thus, a conformal 
type of metric fluctuations can also be considered as the effect of fluctuating mass as in 
the scalar field. For an isotropic, conformal type of stochasticity v{x) defined above, the 
interaction term can be manifestly given as 

^ _ _ 
Cdi = -v{x)iIj{x) yjip{x) = mDv{x)ilj{x)ip{x), (13) 

where we used the equation of motion to obtain the second expression. The equation of 
motion in this model will be 

{i ^ — rriD — vci{x)) ipi^x) = 0, (14) 
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where vd{x) = mDv{x). The corresponding equations for Green functions are 

— (□ + m| + vs{x))G{x, x') = S'^i^x — x') for a scalar field, (15) 
{i ^ — mj:, — vd{x)) S{x, x') = (5'^(x — x') for a Dirac field. (16) 

We expect that the autocorrelation functions of the stochastic fields have the following 
forms: 

As{x -y) = {vs{x)vs{y)) = usJ^lUx - x'\/L)J^T{\t - t'\/T), 

Aoix-y) = {vDix)vDiy))=UDJ'Li\x-x'\/L)J^Ti\t-t'\/T), (17) 

where J-'l and J^t are rapidly decaying functions and L and T are characteristic correlation 
ranges of the stochastic field in space and time, respectively. L and T can be regarded as 
minimum length and time appeared as a result of coarse graining of microscopic dynamics 
characterizing the finite resolution of spacetime. The functions JF^ and J-'t are normalized 
such that / J^id^x = J J-'xdt = 1. For a possible choice of the form of fluctuations, 

{v{x)v{y)) = J"'^^ e-l"-"'l/'^' (18) 
\x — x'l"" 

for a constant a, where L = Ipi = 1.6 x 10~^^ cm is the Planck length, we have 

Us = 47rm^/J, (19) 

independent of a. The disorder-averaged retarded (advanced) Green function {Gr(^a){p,p')) = 
Gr{a) {p) can be written in the form 

Gb.{A){p) = — 2-——- — / \^ / \ ^ (20) 

where mn is a renormalized mass and S/(p) is an imaginary part of the self energy. The 
upper (lower) sign corresponds to the retarded (advanced) Green function. For the stochastic 
field that obeys Eq. (0), v{x) can be approximated as a white noise potential if |p| << Mpi, 
where Mpi = 10^^ GeV is a Planck mass. Consequently, the self energy S/(p) is independent 
of the momentum and depends only on the frequency in such a limit. That is, the low energy 
elastic scattering is insensitive to the length scale that characterizes the fine structure of the 
medium. In this limit, the effect of the medium is to renormalize the frequency of the wave 
propagating through it. The real part of the self energy can be absorbed in the frequency 
term in this limit and will not be considered in this paper. Such a medium is called the 
effective medium. Under this condition, the mean free path of this system is given by 
Im = 4:71 /us in the Born approximation. Combining with Eq. (|T^), we obtain Im = i^s^lpi^- 
For a mass rris much smaller than the Planck mass { ms << Mpi), this naive calculation 
illustrates that the mean free path is much larger than the Planck length [Im » Ipi)- Similar 
results hold for the spinor field as well. The coherence length Lcohi on the other hand, is 
determined by many factors including the other possible interactions, the time dependence 
in the randomness, and external fields and will not be specified here. 



6 



3 Mesoscopic Effects in Stochastic Minkowski Space- 
time 



In this section, we consider quantization of the noninteracting scalar field that obeys the 
stochastic equation of motion discussed in the last section. The parallel argument for the 
Dirac field is given in Appendix A. The stochastic action for the scalar field has the form 



Svs [0, 0l 



d X 



(21) 



This action is invariant under the global U{1) gauge transformation 

(f){x) e'^(j){x), 



X e 



The corresponding Noether current is 



id,, 



(22) 



(23) 



3.1 Closed time path formalism for the stochastic system 



Imt A 



Re t 



Figure 1: The contour in the closed time path integral formalism. 

The quantization of the stochastic system described in Eq. (1) can be treated in the 
closed time path formalism. Here we assume that the whole system consists of a complex 
scalar quantum field and a classical stochastic field vs- Then the density matrix for the 
whole system is given by 

P5[0i,0i,02,0Lws,t] = {(f)u (j)\,vs\ps{t)\(l)2, (pLvs) , (24) 

where |0, 0^, vs) is an eigenstate of the field operator for a particular realization of vs such 
that 

^ixM^\vs) = <PixM<p\vs). (25) 
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To obtain the reduced density matrix ps[(pi, 0i, 02, A for the system, we average over vs 
as 



Ps[(pl, 01, 02, 02, t] = {ps[(pl, 01, 02, 02, Vs, t])v 

The closed time path partition function for this system is given by 
Z[J,j\vs] - 



(26) 



(i0/ci0j-(O_|Texp -i J2 • \(j)f,(j)^f,vs) 
X (0/,0j,V5|Texp [iJi ■ 0] |0_) 

= J d(t>fd(t>\D(t>^D(t>\D(t>2D(t)l exp [i(,S[0i, 0l] - ,5[02, 0|] + Ji • 0i - J2 • 02) 



X exp 



i(S'/[0i, 0l, - S'/[02, 02, vs\) 



(27) 



ion psiU) = |0-)(0_| and 5j[0,0t,-i;5] = - j d^x vs(l)^ J-4) = ^d^x [j^{x)(t){x) + (t^{x)J{x) 

and T denotes an anti-time-ordered product. The path integral above is defined along the 
two paths, one forward in time and the other backward in time (Figure 1). For an arbitrary 
initial state given by ps[0ii, 0ii, 02i, 02i], the partition function takes the form 



Z[J, J\ vs] = I d(t>fd(t>)D(t>^D(t>\D(t>2D(t>\ exp [i(5[0i, 0l] - ^[02, 0^] + Ji • 0i - J2 • 02) 

X P5[01i, 0L, 02i, 0L]^[01, 01, 02, 02, Vs]-, 



where 



F[0i, 0|, 02, 01, Vs] = exp \i{Si[(f)i, 0I, Vs] - -S/[02, 0L Vs]) 



(28) 



(29) 



is a stochastic influence functional. 

In the absence of the stochastic field. 



Z[J, Jt] = j d(t>fd4>\D4>^D4>\D4>2D4>\ exp [i(5[0i, 0l] - 5[02, 0^] + Ji • 0i - J2 • 02) 

X Ps[0ii, 0L, 02i, 0L] = exp [-i J^GJ] , (30) 
where now the Green function G acquired the 2x2 matrix structure as 



G 



Gil G12 
G21 G22 



(31) 



and = (<^i,~<^2)- For an initial thermally equilibrium state with the temperature T = 1/P, 
each component of the matrix Green function in the momentum representation is given by 

Giiip) = -GIM = 0{po)Gf{p) + 9{-Po)GUp) - ^ttz sign(po)nB(p)5(/ - m|), 
G'i2(j9) = -27ii sigia.{po)nB{p)S{p^ - ml), 

G2i{p) = -2msign{po)e^^P'^-^'^nB{p)6{p^-ml), (32) 
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where Gf{p) = {p'^ — rrig + ie)'^ is the vacuum Feynman propagator, nsip) = (e^^^^^^^^ ~^)^^ 
is the Bose distribution function, and /i is the chemical potentiaL G can be diagonahzed by 
multiplying matrices ub from each side as G = UBGdU]^ri^ where 



Gd 



Gr 
Ga 



(33) 



where Gr and Ga are retarded and advanced Green functions, respectively. Here ub{p) is 
the thermal Bogoliubov matrix and r] is the 2x2 Lorentz matrix which have the following 



forms [24, 25 



I g-/3(po-M) 
1 



and rj 



1 
-1 



(34) 



Defining the thermal doublet 



<t>^G-\ 



^2) and its conjugate 
= (j)^riUBGa^UB^(f) 

= (j)UBG'^^U]^^(j). 



)J, 02)^7) can write 



By redefining the field variables by global Bogoliubov transformations u^^c, 
(j), one can write the partition function in Eq. (30) as 



(35) 

and (f)UB ^ 



Z[J,J] 



£)0D0exp iSo[(f),(j)] + iJ(f) + i(t)J ps[0j,0i] 



where 



(36) 



Sn 



dPxduj 



(x, uj) {uj'^ + d"^ — m| + ie{uj)i]j 0(x, uj) 



(37) 



Here e{u!) = e x sign(ci;) for some infinitesimal constant e. When the stochastic field vs has 
no time dependence as in Eq. (^, the close time path action is given by 



Syg[(j), (j)] = I d^xduj (j){x, uj) (uj"^ + — m| — Vs{x) + ie{uj)ri^ (j){x, uj) 



Next we express the Green function and the partition function for the scalar field in 
terms of the nonlinear sigma model. After averaging out the partition function in Eq. (28) 
without the source term with respect to the stochastic potential f 5 which obeys the following 
Gaussian probability distribution 



P[vs] = ATexp y" d'^xd'^y vs{x)As^{x - y) vs{y) 



(3^ 
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with the normahzation constant A/", we obtain the reduced action 

Z = {Z[vs\)vs = J DvsP[vs]D(j)D(f)exp iS^g[(j),(f)] 
= J D(j)D(j)exp iSo[(j),(f)\ + iSj[(j),(l)] 

where 

Si[(l>, = ^ / d^xdSHx)(j){x)Asix - y)(f){y)(j){y). 



(39) 



(40) 



Note that for a local correlation As{x — y) ~ 5^(x — y), we obtain the effective 0^ theory 
similar to the one obtained from spacetime foam 0. 

Now we extract the slow modes by the Hubbard- Stratonovich transformation. Introduc- 
ing the auxiliary bilocal matrix field a{x, y) as p9| 



(41) 



1 



Da exp [-- / (iV(i^?/Tr[(T(x,y)A5^(x - y)cr(y,x)]] exp iSusWA,' 



where 



SHs[<y-, I 



d xd y (f){x)a{x,y)(j){y) 



and the trace is taken over thermal indices. The partition function can be written as 

1 



DaD(j)D(f)exp 



Tr[aA^V; 



exp 



(42) 



(43) 



where a is Hermitian by construction, i.e. cr^(a; — y) = (T{y — x). In energy representation, 
Eq. (42) becomes 



ShsW, ' 



n o .dujduj' J,-, N -'•^ /N 

d xd X (p{x, ct;)(Ti^^/(x, x'j(p{x', uo ) 

Itx 27r 



(44) 



and 



5*0 [0,0] + SHs\<y: 

[x 



(45) 



d^xd'^x'^-^ (f){x, to) (uj"^ + — m| + ie{uj)T]] 5{x — x')5{uj — to') — cru;uj'{x, x') 
27r 27r LV / 



{x> Uj') 



After integrating out and 0, we obtain 



Z = 
X exp 



J Da exp 



IjTrlaAs'a] 

,o ,o ,duj duj' , 
d'^xd'^x Tr log 



(46) 



[uj"^ + — m| + ie{uj)i]j S{x — x')S{uj — uj') — a^^i{x, x') 
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The nonlinear sigma model is commonly used in many different areas in physics to study 
collective excitations and the dynamical symmetry breaking property of the system. Never- 
theless it has not been studied previously in the closed time path method in detail. When 
applied to the disordered systems, it is required that all the fermion loop diagrams will cancell 
in order to include the effects of elastic scattering due to impurity which carries no energy. 
Such techniques as replica formalism |[T^ or supersymmetric extension are commonly 

, m. M M 



used for this purpose. A general class of real time path ordered methods ^7 



is also known to have such a property owing to the energy integral. The closed time path 
formalism employed here has the advantage compared to other methods in that it is natu- 
rally extensible to the nonequilibrium setting. The trace part in Eq. (^) can be expanded 
in terms of the a field as 



Tr log {[uj +d — nig + ie{uj)r]] 6{uj — uj')6{x — x') — cXi^uj'iXjX')} 



Tr log 



a 



Trlog \G-^] 



n=l 



n 



Tr[Ga]", 



(47) 



where G{x,x') = (□ + ^«|)"^ is the free boson propagator. 

Now we assume that the spatial correlation of disorder decays sufficiently fast so that 
the a field can be treated as a local field variable in space. This implies that the stochastic 
potential vs has the time- independent form 



P[vs\ = Afexp 



2m. 



d^x Vs{x) 



In this case the partition function is reduced to 

1 



Z = [ Da exr> / 



X exp 



- / d^x 



2us 
duo duj' 
27r 27^ 



d^xTia'^ix] 



(4J 



(49) 



Trlogl^cj^ + (9^ - m| + ie{uj)ri^ 6{uj - to') - (t^^>{x)} 



The equation of motion obtained from Eq. 
coherent potential approximation: 



9D corresponds to the one obtained from the 



Us 



(50) 



The real part of the a field gives the mass and frequency renormalization and will not be 
considered hereafter. The imaginary part of the a field (TiIuj) gives the scattering rate. 
Writing the homogeneous solution of Eq. (^) as <7^^/{P) = —iaj{u!)sign{u!)6{P)6{u — u'), 
one obtains the relation 



usttN^u) 
2M 



(51) 
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where N{uj) is the density of states of the scalar field. 
The quadratic term in Eq. (^) can be written as 



- / d^xi(fyid^X2d^y2 



duj duj' 
27r 27 



Tr 



ba 



In the momentum representation, the above expression becomes 



d'^P d^k 



rTr 



(27r)4 (277)4 

The Fourier component of the a field is defined as 

d^P d^k 



G'\k + ^)al\~P)G'^\k - ^)al\P) 



(53) 



(yk{P)e 



-ikr —iPX 



(54) 



where r 



X 



(27r)4(27r)4~ 

y and X = {x + y)/2. Thus we obtain the kinetic term of the a field in Eq. 



m 



d^p d'^k 

(27r)4 {2tiY 



Ti{al\-P) 



--G''\k-^)G'\k + ^) 



(55) 



If the spatial dependence in the a field is local, the a field is independent of the momentum 
A;, i.e. <Jk{P) = crko{,P) we can integrate out k and obtain 



d^P dh 



(27r)4 27r 



°Tr{<(-P) 



1 

Us 



(G^^k - -)G^\k 

Jk 2 



p. 

2' 



a 



ko 



(56) 



By expanding inside the bracket in Eq. ( [56| ) with respect to Pq and P, one can show that, 
in the limit (Ji{ko) » koPo, the off diagonal term of the free propagator of the a field in the 
thermal indices gives the massless excitation which has the form of the diffusion propagator 
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KiP)<i-p)) 



ba I 



2aj{ko 



7iN{ko) D{ko)P^ - iPoV" 



(57) 



where D{ko) is the diffusion constant. Diagramatically this form is obtained by including 
all the ladder diagrams in the particle-hole propagator. The diffusion constant is related to 
the dc conductivity Co through the Einstein relation: Gq = N{kQ)D{ko). The diagonal term 
in Eq. (5B) gives the massive excitation that can be integrated out. However, it does not 
contribute directly to the infrared divergence responsible for the universal behaviour of the 
conductance fluctuations and will not be considered in this work. 

In this U{2) nonlinear sigma model, the matrix field a takes its value in the coset space 
U{2)/U{1) X f/(l). After making a transformation a V'^aV with V = 9{uj)l + 0{-uj)a2, 
where matrices 1 and a2 act on the thermal indices, the saddle point solution changes its 
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form as sign{uj)ri rj. In this representation, the field a around the saddle point can be 
parametrized as ||2l[| 



a 




-q^q 



(58) 



The matrix fields q = g^^^/ and = ql^,^ carry two frequency indices. Inserting this expression 
in Eq. (^61), we obtain the form of the free propagator of q field as 

{qkoiP)qli-P)) 



(59) 

(Ji{uj) in the numerator was absorbed in the redefinition of q. Thus the partition function 
for the linear part of the sigma model with respect to q field is 

Z = J Dq^Dqexp{-n J J^,^N{ko)ql{-P) [D{ko)P' - ^Po\ g.o(^)}- (60) 

Higher order vertex terms will follow corresponding to the expansion in Eq. (^). The 
fields q and are free from the constraint and take all possible values. Time dependence 
in the stochasticity can be ascribed to the intrinsic property of the effective medium and 
be treated as the frequency dependent random potential. The modification associated with 
this change can be absorbed in the diffusion constant as the term proportional to dvs/duo. 
If this correction is too large, a different approach is required. The usual prescription to 
account for the effect of inelastic scattering is to include inelastic scattering rate A in the 
denominator of the diffusion propagator so that we replace D{ko)P'^ — iPo in Eq. (|60D simply 
by D{kQ)P'^ — iPq + A. A will set the time scale beyond which the phase memory of a 
scattered wave is lost and the transport behavior becomes classical. 



3.2 Mesoscopic fluctuations of scalar fields 

The conductivity associated with the Noether current in the presence of the external field 
with a frequency n can be written by the current-current correlation function by the Kubo 
formula |Q as 

Ck(x, y) = - dujQ^{uj)Y,jmn{x)jnm{y)S{uJ + K- UJn)S{uJ - UJ^), (61) 

K J—oo ^„ 

where jmn = i<Pln d (pn is the Noether current expressed by two energy eigenstates and ^^i^) 
is a smearing function that depends on the characteristics of the system and the environment. 
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Figure 2: The Feynman diagram for the conductance. The metric fluctuations are included 
as a ladder contribution which is represented by the dotted line. 



Near equilibrium, it can be written as ^^^(co') = ps[uj] — ps[uj + n], where ps[uj] is the initial 
density matrix. If the metric fluctuations are independent of temperature as we assume, the 
effect of temperature on the conductivity only arises from this term. Here we assume that 
$1^(0;) is normalized such that / dujQf,{u!) = k. In Figure 2, the conductance is represented 
by the Feynman diagram. In the leading order weak disorder expansion, averaging over 
disorder is taken into account by including all the ladder diagrams. The conductivity can 



be written in the more familiar form in terms of Green functions as p3 



(62) 



1 



duj ^^.iuj) [Gr{x, y, uj) - Ga{x, y, u)] dxdy [Gniy, x,u + k) - GAiy, x,uj + k)] . 



Ann J-c 

With the expression of Green functions in terms of thermal fields, 

GRix,y) = -ie{xo - yo){y){x),(i)\y)\) = i{(f)i{x)(f){{y)), 
GA{x,y) = ie{yo-xo){y){x),(i)^{y)]) = -i{(j)2{x)(l)l{y)), 
we write Eq. as 

Cf,ix,y) = 



(63) 



-1 
Attk 



duj fi«(cu) ^(0"(x, u) dx (p^\x, UJ + K)(f)'{y, uj + K)dy (p^'^iy, uj)). (64) 



abed 



Note that even though the thermal indices in the sum run over all possible values, only 
pairwise equal terms contribute to the conductivity. 

This expression can be obtained directly from the partition function in terms of the 
functional derivative by introducing the external source term A'^ in the form 

duj duj' 



S[A, 



d'^x- 



27r 2n 



(f){x,uj)A'^{x)- d 6^(j){x,uj'), 



(65) 
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with 



= 6{uj' — uj + k) 



1 1 
-1 -1 



(66) 



The matrix in 6^ acts on the thermal indices. Note that 5^ is nilpotent, i.e. 5^ = 0. The 
nonlocal conductivity is given by 



-n 6^W[A] 



(67) 



where V is the spatial volume of the system and the integral over energy indices are under- 
stood. In the presence of the source term, the partition function for the matrix field a in 
Eq. (^61) becomes 



Z[A] = e*^[^] 



Z^crexp 



Tr 



exp 



duj doj' 
27r 27 



(68) 



Tr log{ 



(cu^ + 9^ - m| + ie{uj)7]) 5{uj - uj') + iA^- d 5« 



6{x - x') - (r^^>{x, x')} 



Making use of the gauge symmetry in Eq. (|68D , for a constant vector field A'^, the source 
term in above can be generated by the following gauge transformation 



[X) 

ix) 



[x)e (69) 
Correspondingly, the Green function G(x, uS) and the a field transform as 

G{x,y) ^U~\x)G{x,y)\J{y) (70) 

and 

a{x,y) ^U-\x)a{x,y)U{y), (71) 

where \J{x) = e'^^'^"^"-. This gauge symmetry will induce the gauge coupling in the effective 
Lagrangian through the covariant derivative V = d + iA'^S^: 



(72) 



Tr [yaVaj 

= Tr ((9(T(9(t) + 2«Tr (I'^^a^a) - Tr (^A^' 6 A""^ 6 ^,a) . 
Now we obtain the expression of the conductivity in terms of the a field: 



-TT 



6^W[A] 



— 71 



(73) 



d'^Xi(TT[6i^a{xi)6_i^a{xi)]) - / d^Xid^X2{TT 6^a{xi)da{xi) Tr 6_^a{x2)da{x2) ) 
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and its fluctuation: 



Cl 



TT" 



5^W[A] 



1 
V 



■^(1) _ ^(2) ^ ^(3) 



(74) 



where 
and 



27r2 



K1K2V 
7r2 



Tr 



5«;,Cr(xi)(9o-(fi 



(i a;2(Tr [5KjO-(xi)(5_K,cr(fi)] Tr [5«;2cr(f2)5-K20'(^2)]), 

(i^a;i(i^X2(i^X3(Tr 5«;-^(t(xi)9(7(xi) Tr S^^^(j{x2)d(j{x2) Tr [6^2'^{x3)6_^2a{x3)]), 
d^Xid?X2d?X3d?Xi 



Tr 



^-K^(y{x2)do{x2) Tr 5«,2Cr(f3)(9cr(x3) Tr 5_«2Cr(f4)5f^(54) )• 



Here we further included the energy indices in the definition of the trace as Tr {O) = 
TT / dko/{2'K)'^N{ko)D{ko) J cLPq J2a ^""(^o, Po)- We are interested in the dc conductivity eval- 
uated in the limit ^1,^2 0. Note that Q^i^) is generally a peak function which describes 
a wave packet peaked around the specific mode u = ujq. Furthermore, if we assume for 
simplicity that it is given by the step function: fl^i^) = 1 (for Uq < u < uq + n) and 
(otherwise), then taking the dc limit k — extracts the particular mode Uq. The effect of 
finite temperature T can be viewed as an additional smearing due to the width of Qi^{uj). 
From Appendix B, in the dc limit, the formula above gives the simple result 



^2 



lim Cl 



VMjR 



(75) 



where c = 7.295 ■ ■ ■ is a constant and M/^ is the infrared cutoff of the momentum integral. In 
Figure 3, the corresponding Feynman diagrams for the conductance fluctuations are given. 
Other diagrams that contain crossed diagrams contribute as higher order terms in the weak 
disorder expansion. Here we assume that the region with the fluctuating metric is restricted 
in the finite cube with the edge length L and that the rest of the spacetime is flat. This 
enables us to handle the problem as a scattering process. Taking L as the parameter that lies 
in the mesoscopic scale, the analogy with the electric circuit becomes elucidated. If we use the 
relation Cq = gL'^~'^ for the conductivity in d dimension for the dimensionless conductance 
g, we obtain the conductance fluctuations in terms of the conductivity fluctuations as = 
Q2j^2d-A^ Then from Eq. (^), in three dimension, using Mm ~ tt/L and we obtain 



^2 ^ 



2.322 



(76) 



This value is universal in the sense that it is independent of the amount of stochasticity ini- 
tially assumed and the size of the stochastic region MB, p6[. The conductance is also directly 
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related to the transmission matrix. Indeed one can show that the conductance measures the 
intensity of wave transmission and the fluctuations of conductance correspond to the 
fluctuations of wave intensity. 




(c) 



Figure 3: The conductance fluctuations are represented by Feynman diagrams. The shaded 
regions are diffusion propagators. (a),(b), and (c) contain two, three, and four diffusions and 
correspond to C^-^\ C^'^\ and C^^^ in Eq. (|T^), respectively. 



4 Discussion 

In this paper, we showed the analogy between the field propagation in Minkowski spacetime 
with a small stochasticity in the metric and the wave in disordered systems. While the 
electromagnetic field propagation in a stochastic spacetime is similar to that in a random 
media, the scalar and spinor field propagation was shown to be similar to the electron in a 
disordered potential. Both cases can be treated similarly, however, the following difference 
should be noted. In the former, the field remains massless and the randomness affects the 
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refraction property of light differently depending on the frequency of the wave. In particular, 
low energy scattering is suppressed; In the latter, a random mass causes scattering with any 
energy. Mesoscopic fluctuations associated with wave propagation were characterized by 
the nonlinear sigma model in the closed time path method. We introduced the collective 
fields by the Hubbard-Storatonovich transformation and integrated out the fundamental 
field variables and obtain the nonlinear sigma model written in terms of the collective fields 
only. The conductivity and its fluctuations were expressed by these fields. For the time 
independent or slowly dependent stochasticity, the fluctuations of the dc conductivity were 
shown to be universal and of order unity. The origin of this universality is traced back to 
the infrared divergence due to the Nambu-Goldstone boson which appears as a result of 
symmetry breaking. 

Although the induced effects on the propagation of waves in the presence of metric fluctu- 
ations are themselves of theoretical and observational importance as long as the backreaction 
of the matter fluctuations is small, a self-consistent treatment is necessary for Planck scale 
processes[^]. This line of consideration is important particularly of the metric fluctuations 
produced in the cosmological processes. 

Relativistic quantum field theoretic calculation of the transport coefficients has been 



developed during the past decade |^, ^ The electrical conductivity in the early universe 
controlls the generation of primordial magnetic field which is believed to be the origin of the 
strong magnetic fields presently observed in spiral galaxies |^ , The method developed in 
Sec. 3 also gives the field theoretic basis for the study of, for example, mesoscopic fluctuations 
due to random magnetic fields. 

The coarse graining of microscopic degree of freedom necessarily induces the nonlocal 
correlation in the stochastic fields. Moreover, unitarity in the whole system guarantees the 
relation between the dissipation kernel and the noise kernel in the form of the fluctuation- 
dissipation theorem upon coarse graining. In the present work, the nonlocal, noncommu- 
tative origin of the stochastic fields and the effect of dissipation are ignored and only the 
classical aspects are considered. Possible manifestation of the quantum nature of under- 
lying microscopic gravitational dynamics in the mesoscopic effects remains to be clarified 
||43|| . Fluctuating metric is also relevant to infer possible decoherence effects in the quan- 
tum interference of propagating particles [Q, |4^, |4^. The closed time path method gives a 
suitable framework to discuss such effects. Our results of conductance fluctuations assume 
that the time scale of metric fluctuations is relatively long. In such a case, the time scale 
of fluctuations appears in the coherence time scale to restrict the validity of the arguments 
and diagrammatic calculations based on the coherence between different modes beyond this 
time scale. This prescription is quite successful in explaining many mesoscopic experimental 
results such as electron scattering in a helium gas. Thus, the heavy defects randomly created 
in the phase transition in early universe can be the origin of such fluctuations of the metric. 

Since, we have not specified the origin of stochasticity in the metric in this work, explicit 
derivation of such stochasticity from the fundamental model of gravity is desired. Space- 
time uncertainty proposed in the context of string theory may have similar effects as 
discussed in this work on low energy physics. Branes or other solitonic objects that appear 



18 



in string theory acquire heavy mass in the weak string couphng hmit and become another 
possible source of stochasticity. Along with the possible decoherence associated with fluctu- 
ating metric, mesoscopic effects treated here may have an observable consequence on future 
experiments ||48| , These microscopic origins of metric fluctuations are intrinsically be- 
yond the validity of semiclassical gravity. Therefore the metric fluctuations introduced as 
the modiflcation of semiclassical Einstein equation in this paper possibly capture the essen- 
tial effects of near Planck scale physics on the sub-Planckian scale physics effectively while 
the self-consistency based only on the conventional semiclassical Einstein equation may not 
have a predictive power on such a phenomenon. Clarifying the difference between these ap- 
proaches needs more careful study. Effects of metric fluctuations in cosmological and black 
hole spacetimes are considered by many authors, for example, in [^0|, ^ ^ ^ ^ ^ 



To identify how the mesoscopic effects discussed in this paper manifest themselves in such 
curved spacetimes is of particular interest. These directions are currently in progress. 
Acknowledgement The author thanks Prof. B. L. Hu for various discussions which mo- 
tivated the present work. Prof. V. Frolov for useful comments, and Prof. D. Page for 
explaining his relevant work on the spacetime foam. 



A Mesoscopic fluctuations of Dirac flelds 

In this appendix we consider the quantization of the Dirac fleld which obeys the equation of 
motion in Eq. (0). The action for this system has the following form, 

(Al) 



S^^y[ip,tp] = J d'^x^ix) (i ^ -mo - VDix))^{x). 
This action is invariant under the global U{1) gauge transformation 



ipix) 

Tpi^x) 



ij{x)e-''^, 



and the corresponding Noether current is 



(A2) 



(A3) 



In the absence of the stochastic fleld, the closed time path partition function has the 
form corresponding to Eq. (30) 



Z[J,jt] 



diljfdiljfDilJiDilJiDilj2Dip2exp {{Slipi^i/ji] - S['ijj2,'ip2] + Ji ■ ipi - J2 ■ '^^2) 



X PD[i^ii,i^ii,i^2i,'ip2i] = exp 



-iJ^SJ 



(A4) 



where J ■ ip = J d'^x J'^ {x)'iIj{x) + ■ip^{x)J{x 

S = 



Su S12 
S21 S22 



The matrix Green function S has the form 

(A5) 
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whose components are 



^11 (p) 

Suip) 



^(Po)5'f(p) + d{-Po)S*pip) - 2ni sign(po)nF(p)5(p^ - ml,), 
2711 sign(po)^F(p)(^ + mD)5{p^ — m 



D)^ 



2m sign(po)e^^^° ^''nF{p){]^ + mD)5{p^ - m^), 



(A6) 



where Sp{p) = i^/p — mu + ie)~^ is the fermion vacuum Feynman propagator, np{p) = 
(^g/3(po-M) _|_x)-i is the Fermi distribution function, and J"!" = — j])- ^ can be diagonahzed 
by multiplying matrices up from both sides as 5* = upSdUp^r], with 



Sr 
Sa 



(A7) 



where Sr and jS'a are retarded and advanced Dirac propagators. Here up{p) is the thermal 
Bogoliubov matrix which has the following form 



UFip) = ^Jnp{p)ef'^°- 



I _g-/3(po-/i) 
1 1 



(A8) 



Using the above property, one can write 



^t^-i^ = ^^upSfup^^ 
= 'ijjupSj^Up^'ijj, 



(A9) 



where = {ipi, ijj2) is the thermal fermion doublet and we define its conjugate ip = {ipl, ip2)V- 
By changing the field variables by global Bogoliubov transformations Up^ip — > ip and ipup 
■0, one can write the partition function in Eq. (A4) without the source term as 



Z ^ j dipfdipfDipDipexp i5'o['0,'0] 



(AlO) 



where 



In the presence of the random variable vp, which is assumed to obey the probability 
distribution given in Eq. (38), we average the partition function over vp, and obtain the 
reduced action 



{Z[vd]) = J DvDP[vD]Di>Dijexp[iSy^[iP,i;] 



(A12) 
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where 



d xd yip{x)^lj{x) ADix - y)i/j{y)ip{y). 



By introducing the auxihary bilocal matrix field ao^x^y) as 



(A13) 



(A14) 



Dod exp [- j d'^xd^yTx[aD{x,y)Aj^{x -y)aD{y,x)]]eicp iSusl'^D^^^'ip] 



where 



d xd y i){x)aD{x,y)i){y), 



(A15) 



the partition function can be written as 

■1 



DaoDipDip exp 



TrfcT/jA^V/)] 



exp 



iSo[ij,'ip]+iSHsWD,'ip,'ip] ■ (A16) 



In energy representation, Eqs. ( |A15| ) and ( |A16| ) have the form 



SHs[o-D,^,ij] = - I d^xd^x'^^ i){x,uj)aDuu'{,x,x')^{x' ,uj') 



(A17) 



and 



doj doj' 



dxdx' ihix.uj 

After integrating out ■?/' and we obtain 
1 



(A18) 

id ■ ^ — TJiD + ie?77o) 6{x — x')5{uj — lu') — odu)uj'{x, x') il>{x' , uo'). 



Z = 
X exp 



(A19) 



,0 ,0 , duo duo' , 
d^xd^x' Trio: 



2ti 2ti 

and, for the time independent stochastic field as in Eq. 

^ d^xTia']~,{x) 



co'7° — id ■ 7 — m^) + ie?77o) 5{x — x')5{uo — uo') — (Toujui'ix, x') 



Z = J Dao^^p 
X exp 



2ud 

(i^xTr log{ (to'7° - id - ttld + ie777o) 6{uj - uj') - (Jdww'{x)] 



(A20) 



The equation of motion can be obtained from Eq. ( |A2CI|) by functional derivative with respect 
to (Jd- Following the steps from Eq. (|52D to Eq. (|55D, the kinetic term in the ao field in 
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Eq. ( |A19| ) can be given similarly to Eq. ( |55[ ) which enables us to constuct the effective field 
theory in terms of the collective field. 

By the Kubo formula the conductivity can be written in terms of the Green functions as 
in Eq. (H) 

C.{x,y) = (A21) 
dujn^{uj) Tr {[SR{x,y,uj) - SA{x,y,uj)]^[SR{y,x,uj + k) - SAiy,x,uj + k)]^}. 



Attk, Jo 
We obtain 



Attk Jo 



C^{x,y) = 

abed 



This expression can be also obtained directly from the partition function by functional 
derivative after introducing the source term in the form 



i [ d^x^^ ^p{x,uj)A'' -id^i/jix^uj'), 
J Zrr Ztx 

where A'^ is the external source field and 5^ was defined in Eq. (| 
is given similarly by Eq. (0). The gauge transformation 



iIj^x) 
%l){x) 



ip{x)e^ 



(A23) 

Then the conductivity 
(A24) 



as we saw in Eq. (|69D, generates the gauge coupling in the effective theory represented by 



doj duj' 
'27r 27r 



(A25) 



Z = J DaD^W ^2 J "^^Wd^ ^ctd] exp J d'^xd^y 

Tr log { (u;7° - id ■ ^ - mo + «er?7o) ^i^ - +1^" -^6^, 6{x - y) - aocucu' {x,y)} 

This allows us to write the mesoscopic fluctuations in terms of ctd flelds. We can discuss 
universal fluctuations parallel to the scalar flelds in this formalism. For an initially thermal 
equilibrium state, = Pd[^] —pnl^ + 'i], where Pd[^] = nF{uj) is the Fermi distribution, 

the dc limit k, —>■ extracts the Fermi energy up that reminds us of the electron transport 
problem. 



B Conductance fluctuations 



dudu'n^iiJ - K)n^{u;')N\ko)D\ko) J d^P 



K1K2 (27r 



(B26) 
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For K1K2 — > 0, 



y| 5 1 

C(2) ^ ^ . / dudoj'n^iu; - K)n^{uj')N\ko)D%ko) I d?PP^ (B28) 

'+K2 UJ i-P)) 

+ h.c] . 



For 0, 



C(^) ^ ^ = (B29) 



C(3) = ^ / dwduj'n^{u - K)n^{uj')N^{kQ)D^{ko) I (fPp^ (B30) 

[ P)){q.'+.. U-P)ql .'+.,(^))(9^'^(->)?L'(^)) 

(-P)?i'+«2 + ^-c-l + three other terms. 



For K1K2 — > 0, 



r(3) ^ 2^ 72 
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